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ON THE RADIATION IMPEDANCE OF A RECTANGULAR PISTON 


by 

Harold Levine 

Department of Mathematics 
Stanford University 

Single integral representations for the resistive and reactive components 
of the radiation impedance appropriate to a rectangular piston are established, 
thereby enabling a systematic refinement of estimates at both short and long 
wave lengths. Comparisons with previous analyses are made explicit as well as 
extensions and corrections thereto. 

§1. Introduction 

Acoustical literature contains numerous theoretical analyses pertaining 
to the characteristics of planar pistons which encompass both time-periodic 
and transient states of motion. The circular shape has been dealt with in the 
fullest manner, given the most favorable hypothesis from the standpoint of a 
boundary value problem, namely that the piston is set in an infinite rigid 

and the particularly difficult task here, typical of most radiation or 
scattering problems, involves the development of precise details at very short 
wave lengths. Analogous considerations (of a more recent date) bearing on 
the rectangularly shaped piston differ somewhat in matters of derivation and 
presentation, without leading to a genuinely efficient or practicable scheme 
of analytical estimation in the short wave length limit. Thus, if the piston 
sides are of lengths a , b respectively and its uniform normal motion is 
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time-periodic at the (angular) frequency co = kc , Chetaev (1951) obtains (by 
an indirect and complex reduction of multiple integrals) the resistance and 
reactance formulas (omitting normalization factors) 

I 
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specifies the geometrical aspect ratio and 

= Ccx = -I 

•X ’ 


( 2 ) 


at; 

5 ' 


<um.X 


AX 


designates the sine/cosine integral functions. Burnett (1969) recasts the 
formulas (1) through an integration by parts which eliminates the latter func- 
tions, yielding (in a modified notation) 
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where 


Burnett writes: "Despite the deceptively simple appearance of the integrals 
in (3) , this author was unable to further reduce them in terms of known tab— 
ulated functions. They can be expressed in terms of confluent hypergeometric 
functions but the various forms investigated do not correspond to any known 
functions. Since Chetaev enjoyed the assistance of two very able mathematicians, 
Tikhanov and Samarskii, yet left his solution in the form of (1), it might be 
construed that there is no easy reduction of these integrals." In keeping with 
such observations he proceeds to a numerical evaluation of the representations 
(3) on the basis of an accurate quadrature technique (further explained in 
Burnett and Soroka, 1972) and to develop an analytical technique, suitable for 
large 6 , which "becomes prohibitively tedious" in its execution; no specific 
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estlmates of the quantities R and X applicable when 6»1 are furnished, 
however. Stepanishen (1977) deduces the radiation impedance by complex Fourier 
transformation of the (previously known) impulse response function of the 
piston (which is itself a Fourier transform of the impedance); and notes that 
the ensuing representation can be brought into full accord with the one due to 
Burnett. He relates the impedance to particular (parameter-dependent) values 
of a complex integral 

^ (3) 

o CoJ2^ ^ 

and connects the asymptotic behavior of C , as x , with those of R 
and X for ka >> 1 , The given estimate 

,f»l C3., 


(attributed to an elementary application of stationary phase techniques) 
should involve a (lower order) term, namely 0(x“^) , and this necessitates a 
change in the stated versions of R and X at high frequencies (which is made 
explicit in §3) . Stepanishen also obtains a general series expansion for 

in terms of Bessel functions and notes that the corresponding repre- 
sentations for R , X which follow are conveniently approximated at low fre- 
quencies. A like degree of reduction which leads to integral expressions for 
the mutual radiation Impedance of square and rectangular pistons in a rigid in- 
finite baffle has been achieved by Arase (1964), although the nature of theoret- 
ical estimates at short wave lengths is not considered. 
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It is the purpose of this paper to show that the radiation impedance of a 
rectangular piston admits a representation with single integrals which are 
related to cylinder functions and are appropriately constituted for purposes 
of systematic analytical estimation at both high and low frequencies. The 
straightforward manner of arriving at this representation is detailed in §2 
and estimates obtained therefrom are recorded in §3; a separate account is 
intended for the analogous handling of impedance functions pertaining to rec- 
tangular panels or beams, whose prescribed normal velocities have a non-uniform 
modal form. 


§!• Formulation and Transformation of the Impedance Function 

The half-space radiation impedance, Z , of a piston with a coplanar 
rigid baffle, as inferred from the integral of the time— average product of 
the fluctuating pressure and (uniform) normal velocity over its finite surface 
area S , takes the form 





' 


( 6 ) 


where p , c designate the equilibrium density and sound speed, respectively, 
of the adjacent medium and a complex time factor applies. If S 

has a rectangular shape and 







the quadruple integral (6) can be reduced, via the transformations 
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to the double integral 
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that is usefully resolved into a pair spanning the triangular subdomains of 

the rectangle 0<^<a , 0<^<b whose specifications in polar coordinates 
are 
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after affecting the elementary p-integration this yields 
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with only single integrals remaining. 
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In the limit b -*• “ the above representation implies that 
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on invoking the change of variable C = sec 6 . It is evident that the latter 
result describes the radiation resistance for an infinite strip piston of width a 


given by the analogue of (6), namely 
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in terms of a Hankel function or, alternatively, by 
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(10) 


when use is made of the representation 
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with a contour that passes above/below the branch points at 5 = + k , re- 
spectively, so that 

I r 0 / 

= , 

The equivalence of (9) and (10) is affirmed on noting the relation 
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which expresses the consequence of displacing the original integral into the 
upper half of the ^-plane and allowing for the pole singularity at ^ = 0 . 

When a,b , the respective changes of variable ^ = sec 0 and ^ = esc 0 
convert the integrals in (8) to the forms 
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and 
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It is verifiable at once that the representations for Re(Z/pcab) and 
-Im(Z/pcab) inferred from (14) correspond exactly with the Burnett forms (1) 
of the (unnormalized) piston resistance and reactance. Furthermore, the proto- 


type integral in (13), 






e 





duplicates the one which Stepanishen defines and employs in an equivalent im- 
pedance representation, A straightforward derivation of the hitherto available 
formulas which constitute the basis for theoretical and numerical predictions 
relating to the rectangular piston impedance is thus completed; and a different 
manner of handling same, with patent advantages at short wave lengths, now 
merits description. 


§3. Conversion and Estimation of the Impedance Function 
Let the integral 

>^>1 

I > ' I 


which enters into (13) be rewritten in the composite fashion 
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and consider the one-variable function that figures, 
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with the initial value 
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Sequential differentiations with respect to x yield the relations 
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where (x) is the Hankel function; thus, on integrating (20) and taking 

note of the condition dF/dx -►0,x-*'Oo,it follows that 




Another integration supplies the original function, viz. 
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where the versions (21a) and (21b) are particularly apt. 
leads, in conjunction with the relations 
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to the generally valid result 

+ Y-(HrWS,(^;-H;'(7s)S^W}|^ 

involving Hankel, Hq (x) , and Struve, Sq ^^(x) , functions of the zero^^ 
and first orders. 

For large values of x it is expedient to utilize (21b) and write 
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with the evident capability of securing additional terms that contain higher 
reciprocal powers of x . 


Repeated integration by parts yields, as regards the two-variable function 
of (16), 



and the number of explicit terms which involve (higher and higher) reciprocal 
powers of X (y > 1) can be readily found. 

Accordingly, the impedance integral (15) is realized in terms of a single 
variable component, with alternative and explicit cylinder function forms, and 
a two-variable 'incomplete' cylinder function component; the former, which takes 
into account the full singularity of the original integral at its lower limit. 


can be estimated on the basis of (extensive) cylinder function asymptotics 
while a numerical evaluation of the latter is furthered by writing 


r 




= f .1 




4 • • • 





( 







-15- 


CO 

and expressing / ^ e in terms of tabulated sine/cosine integral functions 

y 

When the simplest asymptotic developments of the Hankel function 

Hq (x) , X >> 1 , are employed in (23) and the latter combined with (16) and 
(2A) the result is 
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only the leading term, 0(x ) , appears explicitly in the estimate (5a) pro- 

posed by Stepanishen. Applying (25) to deduce a short wave length character- 
ization of the impedance from (13) it follows that 
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in the same limit. 


(U)’ 


cov 




(28) 


These expressions for the radiation resistance and reactance of a rec- 
tangular piston incorporate two orders of magnitude, namely and k~^ 


beyond those previously given and correct, in the terms of order k 
analogous approximations of Stepanishen; thus, his contributions 
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the 




and 
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in R and X , respectively, are eliminated by taking into account an order 
1/x term missing from (5a) and present in (2A) [cf. Appendix], 
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Appendix 

In a paper published subsequently to that concerning the rectangular 
piston, Stepanishen (1978) cites the theorem: Let V’(C) be N times con- 

tinuously differentiable in a < ^ < 6 . Let 0<X<1,0<1J<1; then, 
as X +00 , 
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Identifications 



- 18 - 


0^=1^ ^ A = |4.= 1 and «fC^; = 

apply, it is found that 



* 0 (^> 

and the expected agreement with ( 24 ) obtains. 




